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We discuss the theoretical status of inclusive radiative decays of charmonium, with a particular 
emphasis on the QCD description of the photon spectrum, where progress has occurred in recent 
years. We also comment on the possible extraction of a s (Mj/^) and on the possibility to gain 
important information on the nature of J/ij) and ip(2S). 
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I. INTRODUCTION 

Inclusive radiative decays of heavy quarkonium systems to light hadrons have been a subject of investigation since 
the early days of QCD 0, It was thought for some time that a reliable extraction of a s was possible from the 
inclusive *ygg decay normalized to the inclusive ggg decay. However, when the experimental data became available for 
7-H ' J/^ Q, it turned out that the photon spectrum, and in particular the upper end-point region of it, appeared to be 
] poorly described by the theory. The situation was slightly better for the T(15) [4j, where, at least, good agreement 
with QCD was found in the central region [Ej| . In fact, the whole photon spectrum for the T(1S) is now well understood 
thanks to a number of theoretical advances which have taken place in recent years (see @ and references therein). 
Here we will mainly translate in a critical way these advances to the case of the J/ip. 

We will stay in the effective theory framework of Non-Relativistic QCD (NRQCD) 0, Q and Potential NRQCD 
(pNRQCD) 

and our terminology will follow that of [111 ]. Let us remind the reader that heavy quarkonium 
systems enjoy the hierarchies of scales m 3> mv 3> mv 2 and m Aqcd, where m is the heavy quark mass, v -C 1 



the relative velocity of the heavy quarks and Aq CD a typical hadronic scale. States fulfilling Aqcd ^ m v 2 are said to 
be in the weak coupling regime (the binding is essentially due to a Coulomb-like potential) whereas states fulfilling 
Aqcd 3> mv 2 are said to be in the strong Coupling regime (the binding is due to a confining potential). States below 

Qh! the open flavor threshold and not too deep are expected to be in the strong coupling regime whereas deep states are 
expected to be in the weak coupling one. States above (or very close to) the open flavor threshold are not expected 

... to be in either regime. 

> 

^ . IT. THE PHOTON SPECTRUM 

a, 

The contributions to the decay width can be split into direct ( dlr ) and fragmentation (f ra 9) 



dV _ aT dlr + dTf ra 9 
dz dz dz 

We will call direct contributions to those in which the observed photon is emitted from the heavy quarks and fragmen- 
tation contributions to those in which it is emitted from the decay products (light quarks). This splitting is correct 
at the order we are working but should be refined at higher orders, z £ [0, 1] is defined as z = 2E 7 /M, M being the 
heavy quarkonium mass and E 1 the photon energy in the heavy quarkonium rest frame. 



A. Direct Contributions 



The starting point is the QCD formula [1 

dT _ M 
dz 167T 



^ = z-^ImT(z) T(z) = -i I fxe-*" (V Q (nS)\T{J^x)J v (0)}\V Q (nS)) iff (2) 
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where J^{x) is the electromagnetic current for heavy quarks in QCD and we have restricted ourselves to 3 <Si states, q is 
the photon momentum, which in the rest frame of the heavy quarkonium is q = (q + , q~,q±) = {zM, 0, 0), q± — q°±q 3 . 
Different approximations to this formula are necessary for the central (z ~ 0.5), lower end-point (z — > 0) and upper 
end-point [z — ► 1) regions. 



1. The central region 



For z away from the lower and upper end-points (0 and 1 respectively), no further scale is introduced beyond those 
inherent of the non-relativistic system. The integration of the scale m in the time ordered product of currents in ((2|) 
leads to local NRQCD operators with matching coefficients which depend on m and z. At leading order one obtains 

l^o_2-, , z(l-z) , 2 l^ ln(1 _^_ 2 Il^l! lll(1 _ z)j (3) 



T dz z (2-z) 2 z 2 v ' (2-z) 

where 



r 32 2 2 (VQjnS^O^S^VQinS)) 

T = -aa g e Q - , (4) 



and eQ is the charge of the heavy quark. The a s correction to this rate was calculated numerically in Ref. [l3| for the 
bottomonium case. A reasonable estimate for charmonium maybe obtained by multiplying it by a s (2m c )/a s (2mt). 
The expression corresponding to (g]) in pNRQCD is obtained at lowest order by just making the substitution 

(V Q (nS)\0 1 ( 3 S 1 )\V Q (nS)) = ^|i?„ (0)| 2 , (5) 

where R n o(0) is the radial wave function at the origin. The final result coincides with the one of the early QCD 
calculations [j], Q • The NLO contribution in the weak coupling regime reads [1] , 



<*Tnlo = c , ( 3 S \ (VQjnS^V^S^VQinS)) 

dz l[ - 11 m 4 [ ' 

ant it is v 2 suppressed with respect to ([3]) . The new matrix element above can be written in terms of the original one 

HI 



(VQinS^V^S^lVQinS)) _ (M~2m-S 1 /m\ (y Q {nS)\0^ S^VqinS)) 



J m 2 ^ \ II 



In the weak coupling regime £i/m is absent [15j, but in the strong coupling regime it must be kept (E\ ~ Aq CD is 
a bound state independent non-perturbative parameter) . The matching coefficient can be extracted from an early 
calculation [l6[ (see also [ijj )• It reads 



c; (-4) - -p«44 ((aw + \f wM ) \ + ±^f) m 

where Fb{z) and Fw{z) are defined in ref. p^ 1 . 

In the weak coupling regime the contributions of color octet operators start at order v 4 . Furthermore, away of the 
upper end-point region, the lowest order color octet contribution identically vanishes [Hj]. Hence there is no l/a s 
enhancement in the central region and we can safely neglect these contributions in this case. However, in the strong 
coupling regime the color octet contributions may become order v 2 and should be kept at NLO. 

Then in the weak coupling regime (if we use the counting a s (m) ~ v 2 , a s (ma s ) ~ v) the complete NLO (v 2 
suppressed) contribution consists of the a s correction t o (f3l) , the relativistic corrections in ([6]) and the corrections to 
the wave function at the origin up to order a 2 (ma s ) (20l.l2l|. Using m = m c — 1.6 GeV , M = 3.1 GeV , a s (2m c ) — 0.23 
and a s (ma s ) = 0.4 we obtain the solid green curve in Fig. [1] 



1 The last term in J5J was missing in [f| , see footnote 4 of |18| . 
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2. The lower end-point region 



For z — > 0, the emitted low energy photon can only produce transitions within the non-relativistic bound state 
without destroying it. Hence the direct low energy photon emission takes place in two steps: (i) the photon is emitted 
(dominantly by dipole electric and magnetic transitions) and (ii) the remaining (off-shell) bound state is annihilated 
into light hadrons. For z very close to zero it has a suppression ~ z 3 with respect to Tq (see 0, [23[ for a recent 
analysis of this region in QED). Hence, at some point the direct photon emission is overtaken by the fragmentation 
contributions QQ — > ggg — > ggqqj In practise this is expected to happen somewhere between 0.2 < z < 0.4, 

namely much before than the z 3 behavior of the very low energy direct photon emission can be observed, and hence 
we shall neglect the latter in the following. 



3. The upper end-point region 

In this region the standard NRQCD factorization is not applicable [l2j]. This is due to the fact that small scales 
induced by the kinematics enter the problem and have an interplay with the bound state dynamics. In order to study 
this region, one has to take into account collinear degrees of freedom in addition to those of NRQCD. This can be 
done using Soft-Collinear Effective Theory (SCET) [25|, [26[ as it has been described in [U [28| . This region has only 
been considered in the weak coupling regime, which we will restrict our discussion to. The color octet contributions 
are only suppressed by v 2 or by 1 — z. Since their matching coefficients are enhanced by l/a s (m), they become as 
important as the color singlet contributions if we count a s (m) ~ v 2 ~ 1 — z. The formula one may use for the 
semi-inclusive width in the end-point region, which was successful for the bottomonium case, reads 



dz 



dT 



cs 



dr 



CO 



dz 



dz 



(9) 



where CS and CO stand for color singlet and color octet contributions respectively. 

For the color singlet contribution one may use the expression with the Sudakov rcsummcd coefficient in rcf. 



1 dF 



cs 



Tn dz 



Q(M-2mz)— 



3(1(1 



1 

/•(") 

'5/2 



r( Mc ) 2 ^V/5o_ 7 («) r(Aic) 2A-V/3o 



+ 



3/ 



(n) 
3/2 



8[/, 



(")l2 
5/2J 



A 2 



2A l _"V/3 



where the definitions for the different functions appearing in (JTUJ) can be found in 
For the color octet contributions we use 



(10) 



dl 



CO 



dz 



(Hu) a s (fj, h ) e\ 



/I6M0A [2 
V 9m 4 ) 



C(x — z)Ss+p{x)dx 



(11) 



fx u ~ mv 2 and (ih ~ Tn are the ultrasoft and hard scales respectively. C(x — z) contains the Sudakov resummations 
of ref. • The (tree level) matching coefficients (up to a global factor) and the various shape functions are encoded 
111 S s +p(x). See [6| for a precise definition of these objects. 

We would like to comment on the validity of the formulas above. This is limited by the perturbative treatment of 
the collinear and ultrasoft gluons. For the collinear gluons, entering in the jet functions, we have lGeV < My/1 — z, 
which for J/ip implies z < 0.9. The formalism is not reliable beyond that point. For the ultrasoft gluons, entering in 
the shape functions (Ss+p(x)), we have lGeV < M(l — z), which implies z < 0.7. Hence, due to the latter, we do not 
really have a reliable QCD description of the upper end-point region for charmonium. However, for the bottomonium 
system, the shape function above turns out to describe very well the data even in the far end-point region, where it is 
not supposed to be reliable either. In view of this, we believe that the formulas above may provide a reasonable model 
for the description of the region 0.7 < z < 0.9. The outcome is the blue dot-dashed curve of Fig. [T] (the flattening of 
the curve for z > 0.85 is an artifact, see footnote 2 of Ref. @) 
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FIG. 1: Direct contributions in the weak coupling regime. The solid green line corresponds to the calculation for the central 
region at NLO, which should be reliable up to z < 0.7 . The blue dot-dashed line corresponds to the calculation for the upper 
end-point region, which is expected to provide a reasonable model for 0.7 < z < 0.9. The red dashed line is the curve obtained 
by merging. 



4- Merging the central and upper end-point regions 

As we have seen, different approximations are necessary in the central and upper end-point regions. It is then not 
obvious how the results for the central and for the upper end-point regions must be combined in order to get a reliable 
description of the whole spectrum. When the results of the central region are used in the upper end-point region, 
one misses certain Sudakov and Coulomb resummations which are necessary because the softer scales Myl — z and 
M(l — z) become relevant. Conversely, when results for the end-point region are used in the central region, one misses 
non-trivial functions of z, which are approximated by their end-point (z ~ 1) behavior. In Q the following merging 
formula was proposed, which works reasonably well for bottomonium, 

1 oT dir 1 dT c ( 1 oT e 1 dT e 
Tq dz Tq dz \To dz Tq dz 

\ c means the expansion of the end-point formulas when z approaches the central region. This expansion must be 
carried out at the same level of accuracy as the one we use for the formulas in the central region. 

Putting all the ingredients together in formula (TT^J) we obtain the red dashed line in Fig. Q]for the direct contributions 
to the photon spectrum. Note that a deep is generated for 0.8 < z < 0.9 which makes the decay width negative. This 
happens in the region 0.7 < z where the calculation of the shape function is not reliable. A deep was also generated 
in the T(15) case, but the effect was not so dramatic there Q. We conclude that, unlike in the T(15) case, the 
limitations in the theoretical description of the end point region make the merging procedure deliver an unsatisfactory 
description of this end point region for J/ip. Clearly, further work is required to understand better the end-point 
region, in particular a description assuming J/ip in the strong coupling regime would be desirable. For the present 
analysis, one should only use the outcome of the merging procedure for z < 0.7, if at all. Indeed, an alternative way 
to proceed would be to ignore the end-point region and only try to describe data in the central region with the QCD 
formulas for this region given above. 




B. Fragmentation contributions 



The fragmentation contributions can be written as 

dr fra g f i dx 



dz 



J2 J ^C a {x)D ai (^,M), (13) 



where C a represents the partonic kernels and D ai represents the fragmentation functions. The partonic kernels can 
again be expanded in powers of v [n| 

c . = E6 , e] <™« (14) 



-0.00001 



FIG. 2: The red dashed line corresponds to the direct contributions (merging), the blue dot-dashed line to the fragmentation 
contributions (v scaling for Os( 3 Si)) and the solid green line is the total. 

where Q stands for NRQCD operators, C [Q] for their matching coefficients, and <Iq for their dimension. The leading 
order term in v is the color singlet rate to produce three gluons (Q = Oi( 3 Si)). The color octet contributions have a 
1 1 a s enhancement. In the weak coupling regime, which we will assume in the following, they are w 4 suppressed, but one 
should keep in mind that in the strong coupling regime they may become order v 2 . Then the color singlet fragmentation 
contribution is of order a 3 D g ^, 1 and the color octet fragmentation are of order v A a 2 D g ^, 1 (Q = Os( 1 5'o), Os{ 3 Pj)) 
or v 4 a 2 D q ^ 7 (Q = Og( 3 5'i)). We can use, as before, the counting v 2 ~ a s to compare the relative importance of the 
different contributions. The existing models for the fragmentation functions [3(| show us that D q ^ 7 is much larger 
than Dg^-y. This causes the v A a 2 D q ^ 1 of the Os{ 3 Si) contribution to dominate in front of the singlet a 3 D g ^ and the 
octet v cx 2 Dg^~[ contributions. Moreover, the ot s corrections to the singlet rate will produce terms of order ct^Dq^^ ^ 
that is of the same order as the octet Os( 3 Si) contribution, which are unknown. This results in a large theoretical 
uncertainty in the fragmentation contributions, which would be greatly reduced if the leading order calculation of 
C q [Oi( 3 Si)] (this requires a tree level four body decay calculation plus a three body phase space integral) was known. 

For the quark fragmentation function we will use the LEP measurement [3lj and for the gluon fragmentation 
function the model [321. These are the same choices as in (28[. For the Os^So) and 08( 3 -fo) matrix elements we will 
use our estimates in [331 ] 

(J/rl>\O fi ( 1 S )\J/rl>)\ IJ ^ M ~ 0.0012 GeV 3 (15) 
(J/^|O 8 ( 3 P )|J/V>}| M=M ~ 0.0028 GeV 5 (16) 

The estimate for (j/tplOsi 1 So )\J/tp) is compatible with the lattice results [HI (nrqcd and Coulomb algorithms). For 
the 0$( 3 Si) matrix element the same lattice calculation gives 

(j/V>|O 8 ( 3 S0|J/V} nrqcd = 0.0005 GeV 3 < J/^|0 8 ( 3 5r)| J/^) coulomb = 0.0002 GeV 3 (17) 

which is much smaller than using the NRQCD v scaling 

(J/V»|0 8 ( 3 Si)|J/V») ~ v A {J/il>\0i?Si)\J/il>) ~ 0.05 Ge^ 3 (18) 

with v 2 ~ 0.3. With the choices above we obtain the blue dot-dashed curves in Fig. [2] (v scaling) and Fig. [3] (lattice, 
nrqcd algorithm) for the fragmentation contributions. These curves turns out to be very sensitive to the value 
assigned to ( J/ip\Os( 3 S\)\ J/tp). When we put together direct (red dashed curve in Figs. O and [3]) and fragmentation 
contributions we obtain the solid green curves in Fig. [2] and [3j if the merging formula is used, and Fig. [4] and [5l if 
only the central region is taken into account. The shape of this curve in Fig. [3 is in qualitative agreement with the 
early Mark II results for 0.4 < z < 0.7 0. 

III. EXTRACTION OF a B {M J/i} ) 

As mentioned in the introduction, a s {Mji^,) can in principle be extracted from the ratio T(J/tp — > 7direct-^0 over 
r s trong(>//V — * X), X stands for light hadrons, 7dircct for photons produced from the heavy quarks and r str0 ng for 
subtracting from Y the decays mediated by a virtual photon. This maybe done in an analogous way as it has been 
recently carried out for bottomonia in [35| . In order to obtain r(J/-0 — > 7direct-^0 it is important to have a good QCD 
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FIG. 3: The red dashed line corresponds to the direct contributions (merging), the blue dot-dashed line to the fragmentation 
contributions (lattice, nrqcd algorithm, for Os( 3 Si)) and the solid green line is the total. 



FIG. 4: The red dashed line corresponds to the direct contributions (central region), the blue dot-dashed line to the fragmen- 
tation contributions (v scaling for Os( 3 5i)) and the solid green line is the total. 



description of the photon spectrum since one should be able to disentangle fragmentation contributions from direct 
ones. This may be done by restricting the fit to data of the QCD expression for direct contributions to the upper 
end-point region and the part of the central region where fragmentation contributions are negligible. As we have seen 
in the previous section, we do not have at the moment a good QCD description of the upper end point region for 
J/ip and hence a model, like the one in [36j| (see also [37|), might be unavoidable. This expression is then used to 
interpolate data to small z and hence to be able to obtain the full inclusive width for direct photons. Then a s (Mj/^,) 
may be extracted from 



T(J/ip — > 7dircct X) _ 36 e 2 c a 
r s tron g (^/V' —>■ X) 5 a s 



l + 0{a s )+0( 



a s (m) 



) + 



(19) 



FIG. 5: The red dashed line corresponds to the direct contributions (central region), the blue dot-dashed line to the fragmen- 
tation contributions (lattice, nrqcd algorithm, for Og( 3 Si)) and the solid green line is the total. 
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where e 2 = 4/9, the 0(v 2 ) cancel in the ratio, and • • • stand for higher order contributions. In the extraction of a s 
from bottomonium of [35|, 0(a B (m)) corrections were taken into account but not C(^-) which are of the same order 
if a s (m) ~ v 2 and in practise turn out to be very large. These have been included in (38j . 



IV. LEARNING ABOUT THE NATURE OF J/ip AND ip(2S) 



It has recently been shown that if two heavy quarkonium states are in the strong coupling regime then the ratio of 
their total photon spectrum in the central region is predictable from QCD at NLO [39( . If the spectrum of both J/tp 
and tl>(2S) is measured and fits well with the strong coupling regime formula, it would indicate that both J/tp and 
ip(2S) are in the strong coupling regime. Unfortunately if it does not, we will not be able to learn much, because it 
may be due to the fact that J/ip is in the weak coupling regime or to the fact that tp(2S) is too close to the open 
flavor threshold for the strong coupling regime to hold for it (or to both). 



V. CONCLUSIONS 



A new measurement of the inclusive photon spectrum for radiative J/tp decays would be of great interest since it 
has only been measured before by the Mark II collaboration more than 25 years ago. The theoretical progress which 
has occurred since may allow, among other things, for a sensible extraction of a s (Mj/^). The additional measurement 
of the photon spectrum for ip(2S) might shed some light on the nature of these states. No theoretical analysis are 
available for other states like ?y c s or x c s. Experimental measurements would definitively trigger them. 
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